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A note on observational signatures in superluminal unified dark matter models 

Yuko Urakawa* and Tsutomu Kobayashft 
Department of Physics, Waseda University, Okubo 3-4-1, Shinjuku, Tokyo 169-8555, Japan 

We explore the possibility that the dark matter and dark energy are mimicked by a single fluid 
or by a single fc-essence-like scalar field. The so called Chaplygin gas unified dark matter models 
can reproduce the observed matter power spectrum by adding a baryon component. It has been 
argued that the evolution of the baryon fluctuations is particularly favoured for the "superluminal" 
case where the sound speed of the Chaplygin gas exceeds the speed of light at late times, as well 
as for the models with the negligibly small sound speed. In this note we compute the integrated 
Sachs-Wolfe signal in the Chaplygin gas models, focusing on the superluminal case which has not 
been investigated before because of the premature understanding of causality. It is shown that 
the superluminal model leads to large enhancement of the integrated Sachs- Wolfe effect, which is 
inconsistent with the CMB measurements. 

PACS numbers: 98.80. Es, 98.80. Bp 



I. INTRODUCTION 

Current cosmological observations indicate that a large 
fraction of the Universe is made of unknown "dark" com- 
ponents: dark matter and dark energy. Their identities 
and physical properties are still masked in mystery. Is 
dark energy a cosmological constant or an evolving scalar 
field? Can modified gravity be an alternative to dark 
matter and/or dark energy? Are dark matter and dark 
energy really distinct components? In this note, this last 
possibility is explored. Our goal is to reevaluate the via- 
bility of (a particular class of) unified dark matter/dark 
energy models from the point of view of cosmological ob- 
servations. 

Unification of dark matter and dark energy is an at- 
tractive idea because it would be simpler than to con- 
sider two separate dark components. Unified dark mat- 
ter (UDM) models are often described by a fc-essence-like 
scalar field <f> the Lagrangian of which is of the form 
C = p(X,4>) with X := -(dcj)) 2 /2. The scalar field is 
assumed to mimic both dark matter and dark energy [2- 
5]. The basic requirements for a UDM model to be vi- 
able are: (i) the background evolution is composed of the 
early-time matter-dominated phase and late-time accel- 
erated phase; (ii) the evolution of cosmological pertur- 
bations is consistent with observations. Among various 
models we are interested in those which are described 
by [6-8] 
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with X < Vq. In the case of a = 1 this reduces to 
the Dirac-Born-Infcld Lagrangian that appears in string 
theory [9, 10], but it is phenomenologically interesting to 
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consider general a (> 0). The scalar field model (1) has 
an equivalent description in terms of a barotropic fluid 
with the equation of state 



Pch = -Ap, 



ch- 
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The fluid version is referred to as the generalized Chap- 
lygin gas (gCg) model [6-8] . Cosmological consequences 
of the gCg model have been extensively investigated so 
far [11-25]. 

The behavior of cosmological fluctuations basically de- 
pends on the sound speed c s of a fluid, and in the gCg 
case the sound speed crucially depends on the parameter 
a. Thus, the allowed parameter space is strongly con- 
strained by the data on large scale structure (LSS) [26- 
28] to be a < 1, leaving the models which are essen- 
tially identical to the usual cold dark matter model with 
a cosmological constant (ACDM), in particular in the 
simplest modeling of the universe filled only with the 
gCg [13]. The situation seems to be drastically cured 
at least in the linear regime by adding a baryon com- 
ponent [15, 16, 29, 30], as baryon fluctuations grow in 
a standard way to reproduce the observed matter power 
spectrum consistently. 1 Since the linear growth factor 
of the baryonic fluctuation is still suppressed to about 
90% compared to the ACDM model, measurements of 
the cosmological parameter as are expected to yield a 
strong constraint on the gCg model [13]. However, the 
limits from ag are still unclear because as is strongly af- 
fected by non-linear effects [29, 30]. 

Most of the literature so far has focused on the mod- 
els with a < 1 probably because the sound speed ex- 
ceeds the light velocity for a > 1 [31, 32]. However, 
as argued in [33], the superluminal propagation of a k- 



1 While the Chaplygin gas fails to cluster due to its increasing 
sound speed, the baryons can cluster because of the absence 
of the direct interaction between the Chaplygin gas and the 
baryons. This implies that an unusual bias is required to relate 
the density fluctuations of these two components. 
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essence-type scalar field does not necessarily lead to the 
causality violation. Interestingly, the authors of Rcf. [29] 
found recently that the growth of perturbations is partic- 
ularly favoured for a > 3 as well as for very small a (see 
also [34]). In [29], a particular solution was proposed 
to assure the preservation of causality. Their point is 
that the necessary condition to preserve causality is that 
the signal velocity does not exceed the light velocity in 
vacuum [35]. Modifying the Lagrangian (1) in an un- 
observable regime, they showed that the signal velocity 
does not exceed the light velocity. 

This note also focuses on the parameter range a > 1. 
We investigate the integrated Sachs- Wolfe (ISW) signal 
in the case where the sound speed of the gCg can be- 
come superluminal, in addition to the power spectrum of 
the baryon fluctuation studied in [29]. It is known that 
a finite sound speed enhances the ISW effect [4], due to 
which the gCg with 0.01 < a < 1 is excluded [11, 14, 17]. 
We shall show that the superluminal models also yields 
large enhancement of the ISW effect, which clearly incon- 
sistent with observations. In contrast to weak lensing, the 
linear calculation of the ISW effect is reliable, and hence 
this is probably most convenient way of constraining the 
superluminal gCg model. 



II. BASIC EQUATIONS 
A. The background evolution 

We consider the universe filled with the gCg and bary- 
onic components. The two components are assumed 
to interact only gravitationally. As was pointed out 
in [16, 17, 29], it is important to include the baryons to re- 
produce the observed matter power spectrum. Since the 
modification made to cure the causality problem does not 
affect the cosmic evolution until now, we use the equa- 
tion of state (2). It follows from Eq. (2) and the energy 
conservation equation, dpch/di + 3H(pch + Pch) = 0, 
that 
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where t is the proper time, a(t) is the scale factor of the 
universe, and H := dlna/dt is the Hubble parameter. 
The present value of the scale factor ao is set to unity. 
Pch.o is essentially an integration constant and represents 
the present energy density of the gCg, and A:— A/ p c ^. 
The ratio between the pressure and energy density, wch, 
is given by 



dark energy phase, wch — — 1> at late times. The transi- 
tion from the matter dominant phase to the accelerated 
phase proceeds more rapidly for larger a. (The transition 
time is determined through the model parameters a and 
A, as explained shortly.) Note that in the limit a — »■ 
the gCg represents usual dark matter plus a cosmological 
constant. 

Assuming the spatially flat universe, the Friedmann 
equation can be recast in 
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where Ha = 100 hT 1 kms -1 Mpc" 1 is the present Hubble 
parameter and fib. o is the present baryon density frac- 
tion. Let us define the transition time as the time at 
which the deceleration parameter becomes zero. The red- 
shift at the transition time is then given by 
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Instead, to avoid the complication we may use the rela- 
tion 
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to define the transition redshift. Since f2b,o is very 
small, the two definitions make practically no difference. 
Let us then define the present value of the total (non- 
relativistic) matter density, f2 m .o- Noting that 
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for a — » 0, we may read off from this that 

On,o := a,o + (1 - n b ,o)(l - ^) 1/(1+a) . 
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Finally, we introduce the key quantity in the present 
note: the sound speed of the gCg, 
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The sound speed plays a crucial role in the following 
perturbation analysis. The behavior of c 2 is plotted in 
Fig. 1. The sound speed remains small before the tran- 
sition time. For a > 1, c s is particularly small in this 
regime. It starts to grow around the transition, and (? s 
can exceed the speed of light at late times for a > 1. We 
thus call the gCg with a > 1 "superluminal," though c 2 s 
can be larger than 1 only at late times and the model 
with a ~ 1 is in fact "subluminal" because liochl < 1- 
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Equation (4) clearly shows that the gCg indeed mimicks 
both the dark matter and dark energy: we have the mat- 
ter dominant phase, w;ch — 0, at early times and the 



B. Cosmological perturbations in the gCg model 

In the linear analysis, we adopt the longitudinal gauge 
in which the metric is written as 
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FIG. 1: The sound speed as a function of the redshift z for 
z t i — 0.8 and different a. The models with a > 1 show the 
nontrivial behavior: c? s is smaller for larger a in the past, 
which rapidly grows at z ~ 1 to be <? a ~ a at 2 = 0. 



FIG. 2: The evolution of the baryon density perturbation for 
k = 0.01 feMpc -1 . The parameter of the background is fixed 
by choosing z tI = 0.8 and Qb,o = 0.04. The evolution is not 
sensitive to a. 



Here, we neglect the anisotropic stress so that the grav- 
itational potential coincides with the curvature pertur- 
bation in this gauge. We do not consider the entropy 
perturbation, and write the pressure perturbation of the 
gCg as 



Spch = c 2 s Sp C h, 



(12) 



where Spch is the density perturbation of the gCg. 

Combining the (t,t) and components of the per- 
turbed Einstein equations to eliminate Spch, we obtain 
(in the Fourier space) 

$"+ [4 + 3cf + (lnif)'] *' + [3(l + c^) + 2(lnff)'] $ 
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where Sb '■= Spb/pb is the baryon density perturbation 
and the prime denotes the differentiation with respect 
to In a. For brevity, we abbreviate the suffix k if not 
necessary. The energy-momentum conservation of the 
baryonic component implies 



(S h - 3$)' = ^v h , 

(«»b) + -Jj = 0, 
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where Vb is the velocity perturbation of the baryon fluid. 
These two equations are combined to give 
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where we used Eq. (13) to eliminate Equations (13) 
and (16) are sufficient for predicting the ISW effect and 
matter power spectrum in the gCg model. 



III. OBSERVATIONAL SIGNATURES 

A. Evolution of cosmological perturbations 

Let us begin with over-viewing the evolution of the 
metric perturbation $ and the baryon density perturba- 
tion 5b governed by Eqs. (13) and (16). At early times 
we have wch — and hence c 2 s ~ 0. Thus, for z ^> z tr , 
we can find the same solution as in the CDM dominated 



universe. Noting that H 2 o 
and (16) are solved to give 



in this limit, Eqs. (13) 
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where the decaying mode has been ignored. Here, Ck 
is a constant, and a; and Hi are the quantities eval- 
uated at some initial time. In practice, we may im- 
pose the initial condition at the last scattering surface: 
a,; = (1 + zls) 1 ~ 10~ 3 . The primordial power spec- 
trum is encoded in the ^-dependent constant Ck- Having 
specified the initial condition, we separate the evolution 
of perturbations for a > a; from the initial amplitude Ck 
and write the solution to Eqs. (13) and (16) in terms of 
the transfer functions as 

$ fc (a) - C fe T*(fc,a), 5 hk {a) = C k %(k,a). (18) 

The most important ingredient that determines the 
behavior of fluctuations is the magnitude of the sound 
velocity c 2 , or, equivalently, the Jeans length of the 
gCg fluid: \ 2 := c 2 J(aH) 2 . The gCg fluctuations with 
k > Aj 1 are supported by the pressure and oscillate, 
rather than grow. This fact puts a very stringent con- 
straint on the unified dark matter model composed of 
a single gCg fluid, as such oscillations hinder the gCg 
component to reproduce the observed matter power spec- 
trum. One can, however, overcome this early drawback 
by adding a baryon component, because the baryon fluc- 
tuations are not affected much by the late-time oscillation 
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ferent a and z tr = 0.8 [29], compared to the observed one [28]. 



of the other component. In other words, the gravitational 
interaction is not strong enough to eliminate the baryon 
fluctuations built up before the oscillatory regime. The 
typical evolution of the baryon density perturbation is 
shown in Fig. 2. Although the different choices of a lead 
to the different magnitude of the gCg sound speed at late 
times, 5b shows a very similar behavior, irrespective of a. 
One should note here that even in the superluminal case 
with a > 1 the baryons can cluster because the sound 
speed remains small at early times. Interestingly, the 
sound speed at early times is much smaller for a > 0(1) 
than for a < 0(1), and hence the former is much closer 
to the standard ACDM model than the latter, at least at 
early times. 

In Ref. [29] Gorini et al. solved the evolution of 8^, 
and computed the matter power spectrum, focusing in 
particular on the superluminal gCg model. They argue 
that the growth of baryon inhomogeneities is observation- 
ally favoured for a > 3 as well as for sufficiently small 
values of a. For the sake of completeness, we have re- 
peated essentially the same calculation as in [29]. The 
power spectra of the baryon component are shown in 
Fig. 3, which are in good agreement with the observed 
one. To rely on the linear analysis, we discard the data 
with k > 0.3/iMpc^ 1 . 

Now we turn to the evolution of the gravitational po- 
tential $. In contrast to the baryon fluctuations, $ is cru- 
cially affected by the late-time increase in c 2 . As shown 
in Fig. 4, the pressure-supported oscillation appears for 
large values of a, leading to a strong time variation of 
$. This is simply because c 2 is larger for larger a at late 
times. Let us focus on the cases with a — 1 and a = 3 
and compare them in more detail. Since the sound speed 
remains c 2 < (aH/k) 2 for shorter time for a = 1 than 



FIG. 4: The evolution of the metric perturbation $ for k = 
0.01 /iMpc" 1 . The parameter of the background is fixed by 
choosing z tr = 0.8 and Qb.o = 0.04. In contrast to the baryon 
density perturbation, the evolution of $ crucially depends on 
a. 



for a — 3, $ begins to oscillate earlier in the former case. 
In the accelerated phase, z < z tI , the sound speed of the 
a = 3 model grows to become larger than that of the 
a = 1 model, which causes the larger time variation for 
a — 3 than for a — 1 in the oscillatory regime. 

B. ISW effect 

Having thus emphasized the characteristic behavior of 
it is important to study the signatures of observations 
that probe more directly the gravitational potential itself. 
An example is weak gravitational lensing, using which the 
authors of [13] claimed that the allowed parameter space 
of the gCg model (with baryons) is very tiny. However, 
their argument is based on the linear analysis despite the 
fact that nonlinear effects are important. Another nice 
example is the ISW effect [36, 37] . A strong time vari- 
ation of the gravitational potential at z < ztr typically 
gives rise to large enhancement of the ISW signal. This 
can put a stringent constraint on the model parameters, 
and indeed the gCg models with 0.01 < a < 1 are ex- 
cluded [11, 14, 17] (see also [4]). In this note we choose 
to pursue this direction, paying a particular attention to 
the superluminal gCg with a > 1. 

Noting that the gravitational potential coincides with 
the curvature perturbation in the absence of anisotropic 



5 



pressure, the angular power spectrum of the ISW effect 
is, for the case with a spatially flat background, given by 
(see [4, 38-40]) 



Cj sw = 4ttT 2 



where 



rno 
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jl(x) is the spherical Bessel function, and rj is the con- 
formal time. For our purpose, we focus on the late time 
ISW effect. The initial amplitude Ck is obtained from 
the primordial spectrum multiplied by the BBKS trans- 
fer function [411: 



fc 3 |C fc | 2 
2tt 2 
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with 
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(22) 

where q := k/(Qx.oh 2 ) and f^x.o incorporates 
Sugiyama's shape correction [42]: 
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FIG. 5: The power spectra of the ISW effect normalized by 
Tq A^/25, which approximately gives the angular spectrum 
for the SW effect, top: z tr = 0.8, bottom: z tI = 1.0. 



^x,o := ^m,o exp -Qb,o - 



2h ^2^,0 



(23) 



We simply assume that the primordial power spectrum 
is the scale invariant Harrison-Zel'dovich one: A^(fc) = 
constant. 

The Hubble Space Telescope Key Project [44] result 
and the WMAP 5yr data [45] tightly constrain the pa- 
rameters h and fib,o- Since the amplitude of the ISW 
effect varies at most few percents within the parame- 
ter region allowed by observations, we use h = 0.7 and 
^b,o = 0.04. The angular spectra of the ISW effect for 
different a and z% T are shown in Fig. 5. The amplitude is 
normalized by Tq A^/25, which approximately gives the 
amplitude of angular spectrum for the Sachs- Wolfe effect 
Cf w . For the standard ACDM model, a = 0, Cj sw is 
roughly 10% of C ; sw . Therefore the amplitude of the flat 
plateau is essentially determined by Cf w . Contrary, in 
the gCg model with a > 0(1), Cj sw is more than ten 
times larger than Cf w , and hence the amplitude of the 
flat plateau is essentially determined by Cj sw . Reflect- 
ing the oscillating behaviour of the enhancement of 
the ISW effect appears also in the superluminal model 
with a > 1. Since the baryon density perturbation 6b 
does not exhibit the distinctive behaviour in the super- 
luminal gCg model, this enhancement would help us to 
discriminate this model. 



Combining the constraints from the CMB anisotropies 
and the LSS, we now discuss the viability of the superlu- 
minal gCg model. From the above argument, one sees 
that, regarding the ISW effect, a significantly smaller 
value of the primordial amplitude A 2 ^ is favoured in the 
superluminal gCg model compared to the standard value 
in the ACDM model. However, in regards to the struc- 
ture formation, the preferred value of A 2 ^ is larger by an 
order of magnitude in the superluminal gCg model than 
in the ACDM model, since the structure formation pro- 
ceeds less efficiently in the gCg model. When we change 
the primordial amplitude A^, ss 2 to rA^ ss 2 , the angular 
spectrum C; and the matter power spectrum P(k) also 
change into rC\ and rP(k), respectively. Here, we de- 
fine A^ ss 2 as the prefered value to reproduce the mat- 
ter power spectrum, which depends on a and ztr- We 
need r < 0.1 in order for the temperature anisotropy 
to be compatible with observations, but then the pre- 
dicted matter spectrum significantly deviates from the 
SDSS data. As an example, the total angular spectra C; 
for r = 1, 0.1, and 0.05 are shown in Fig. 6. For r = 0.1 
the amplitude of the matter power spectrum comes to be 
inconsistent with observations, but even in this case the 
angular spectrum of the CMB still remains larger than 
the observed one. 
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FIG. 6: The power spectra for a = 3, Ztt = 0.8, and different 
values of A^,, compared to the WMAP 5yr data [45]. r is 
defined by r := (An / A^ ss ) 2 . 



Finally, let us note the dependence on the parameters 
a and Ztr- We have computed the ISW signal for a < 10 
and found that these models are inconsistent with obser- 
vations. The angular power spectrum for the ISW effect 
and the matter power spectrum slightly depend on the 
transition time z tr . Our result is robust to change of z tr 
within the observationally allowed region (which is not 
so small) [21, 46]. Indeed, for 0.4 < z tI < 1.2, the su- 
perluminal gCg model with a < 10 cannot be consistent 
both with the LSS and with the CMB anisotropies at the 
same time. 



IV. CONCLUSIONS 

We have investigated the growth of inhomogeneities 
in superluminal gCg cosmology, in which a single fluid 
component mimics both dark matter and dark energy. 
The superluminal gCg model is interesting observation- 
ally because it can reproduce the observed matter power 
spectrum well and theoretically because it does not in 
fact suffer from the causality problem. We have shown 
that the late-time increase of the sound speed affects the 
evolution of the metric potential $ rather than that of 
the baryon density perturbation <5t>- Since the time de- 
pendent gravitational potential generates secondary tem- 



perature fluctuations through the ISW effect, the obser- 
vation of the CMB angular spectrum at low multipoles 
helps to discriminate the superluminal gCg model. We 
have found that, due to large enhancement of the ISW 
effect, the superluminal gCg model cannot explain the 
SDSS data and WMAP data consistently. Here, we fo- 
cused only on the flat plateau of CMB angular power 
spectrum, combined with the matter power spectrum. 
As is the case in the subluminal gCg model [14, 17], nor- 
malizing Ci to COBE at I — 10, it would be possible 
to give stringent constraints from the acoustic peak of 
CMB, too. 

In this note we have only considered a particular class 
of UDM models. Then, a question is: within the con- 
text of UDM models, how generic is large enhancement 
of the ISW effect? It is now clear that the late-time in- 
crease in c 2 s is crucial for the strong time variation of the 
gravitational potential that leads to the large ISW sig- 
nal. Therefore, to avoid enhancement of the ISW effect, 
c 2 s must be sufficiently small almost all the way from the 
past to the present. This point has already been empha- 
sized in [5] (see also [47] ) . For sufficiently small but finite 
Cj, weak gravitational lensing can be used to discrimi- 
nate among UDM models [43] because small c 2 s pushes 
the Jeans length away to small scales. In order to give a 
reliable prediction about weak lensing, it is important to 
address the nonlinear evolution of perturbations in UDM 
models. This issue would be very interesting and is left 
for further study. 

Let us finally comment on the limiting case a — > oo. 
In this limit, the (scalar field model of) gCg behaves like 
Cuscuton proposed in [48]. While the Cuscuton has the 
infinite sound speed, it does not transport any informa- 
tion. It would be very interesting to explore the conse- 
quence of this exceptional case [49] . 
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